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■r^ij- ' A novel two-sided Bitter decoration technique was recently employed by Yao 

et al. to study the structure of the magnetic vortex array in high-temperature 



superconductors. Here we discuss the analysis of such experiments. We show 
' that two-sided decorations can be used to infer quantitative information about 

o ■ 

J-^ ' the bulk properties of flux arrays, and discuss how a least squares analysis of 

the local density differences can be used to bring the two sides into registry. 

d ' Information about the tilt, compressional and shear moduli of bulk vortex 

configurations can be extracted from these measurements. 
PACS: 70.60-w, 74.60Ec 
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I. INTRODUCTION 



The nature of the ordering of the magnetic flux-hne array in the mixed state of high- 
temperature superconductors is a topic of much current theoretical and experimental in- 
terest. Most direct measurements of the microscopic structure of the flux array have been 
obtained via Bitter decoration experiments at low fields [|I|]- 0. The conventional Bitter 
technique employs small magnetic particles to decorate the tips of individual vortices as 
they emerge at one of the sample surfaces. Since in the high- Tc superconductors the vor- 
tices can wander considerably in the transverse direction as they traverse the sample |^ 
and in addition the intervortex interaction at the surface differs from that in bulk 1^,^, it 
is difficult to unambigously infer bulk properties of flux arrays from conventional one-sided 
decorations that measure two-dimensional correlations of fiux-line tips at the surface. Quan- 
titative information on the three-dimensional structure of the vortex lines as they traverse 
the sample can be obtained via neutron scattering, but only very few such measurements 
have been carried out to-date due to their difficulty and cost. In addition neutron scattering 
is usually feasible only at much higher fields than probed by decorations. Very recently Yao 
et al. used a novel two-sided decoration technique to study vortex structure in single 
crystals of BSCCO. These authors have simultaneously decorated both sides of the sample 
and analyzed how the two-dimensional translational and orientational order of the vortex 
array propagates across its thickness. Such "flux transmission spectroscopy" experiments 
are likely to become an important source of insight about vortex matter in the future. In this 
paper we discuss the analysis of these experiments and show how two-sided decorations can 
be used to infer quantitative information on bulk properties of the flux array. In principle 
all three bulk elastic constants of the flux array, the compressional, tilt and shear moduli, 
can be extracted from these measurements. 

Much of the analysis described below was originally carried out for flux arrays in a liq- 
uid phase 0, [||- [0. This case may be relevant to the decoration experiments that are 
field-cooled below the irreversibility line. Because of long relaxation times, the observed flux 
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patterns do not represent the equilibrium configuration of tlie vortices at tlie low tempera- 
ture where the decoration takes place, but may be better approximated by the equilibrium 
configuration at a higher temperature Tf where the flux flux array falls out of equilibrium. 
The value of Tf is not known, but it is estimated to be very close to the experimentally 
observed irreversibility line. Depending on the field strength, the flux array may be in a 
crystalline or liquid-like state when it drops out of equilibrium. We summarize here results 
for liquid, hexatic and crystalline vortex arrays. We also discuss the extent to which config- 
urations on opposite sides of a sample may be brought into registry by a least square fit of 
the difference in the local vortex density. 

We discuss quenched random disorder here only for point pinning and weak surface 
disorder in the flux liquid. The effects of bulk point pinning at the elevated temperatures of 
the low fields irreversibility line are then weak because the impurity potential in thick samples 
is screened out by thermally induced vortex collisions 0. It should be straightforward, 
however, to extend much of the analysis summarized here to other types of bulk and surface 
pinning in crystalline, hexatic and liquid phases. Strong surface disorder could certainly 
obscure the interpretation of double-sided decorations. If surface pinning is not a factor, 
it woud be especially interesting to consider the effect of correlated disorder, in the form 
of columnar pins, either parallel or splayed |l3l, which pass completely through the 



sample. The key experimental question in this case is whether vortices always track a single 
column as they traverse the sample, or if they hop from column to column. This question 
plays a particularly important role in theories of vortex transport in the presence of splayed 
defects [|13| . 



II. TRANSMISSION OF DENSITY FLUCTUATIONS 

Density fluctuations of flux lines in three dimensions are described by the correlation 
function of a coarse-grained areal density field. 
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n{r^, z) = Y, <^(r± - ri{z)), (1) 

i=l 

where ri{z) is the position of the i-th vortex in the {x,y) plane as it wanders along the z 
(z II H) axis. In a sample of thickness L in the field direction and cross-sectional area A, 
translational correlations between the two opposite surfaces of the sample are described by 
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noS{q±,L) = 5n{q±,L) 6n{-ci±,0) > 



-< (5n(qx, L) >< (5n(-q±, 0) >J , (2) 

where Sn{q±, z) = n{q±, z) — rio^'^qx.o denotes the fluctuation of the in-plane Fourier trans- 
form of the coarse-grained fiux-line density from its equilibrium value Hq = B /(pQ. A factor 
of 77-0 has been extracted in the definition of the structure factor so that S{q±,L) — > 1 as 
g_L ^ oo. The angular brackets denote a thermal average and the overbar the average over 
quenched impurity disorder. The subtracted term on the right hand side of Eq. ^ vanishes 
in the absence of quenched disorder. 

Almost thirty years ago Pearl ^ showed that the interaction between the tips of straight 
flux lines at a superconductor- vacuum interface decays as l/r_L at large distances, with r± the 
distance between flux tips along the interface. In contrast, the interaction between fiux-line 
elements in bulk decays exponentially at large distances. For this reason Huse argued that 
at low fields, where the intervortex separation is large compared to the penetration length, 
surface effects may play the dominant role in determining the magnetic flux patterns seen at 
the surface P| . The question of the interplay between bulk and surface forces in determining 



the vortex structure at the surface was addressed by us [|1T| with a hydrodynamic model 
that incorporates the boundary condition on the fiux-lines at the superconductor-vacuum 
interface - which is responsible for the l/r± interaction - as a surface contribution to the 
free energy of the flux array, coupled to the usual bulk free energy. This model neglects all 
spatial inhomogeneities in the z direction other than the presence of the sample boundaries. 
The thermal contribution to the structure factor defined in Eq. ^ was found to be 

STiq±,L) = S2T{q±)Riq±,L), (3) 



where S'2t('?±) is the two-dimensional structure factor at one of the two surfaces and R{q±, L) 
measures the "transfer" of information about density fluctuations across the thickness of the 



sample. For the hydrodynamic model considered in Ref. |Tl| this transfer function is given 
by 

' 52(g±)sinh(L/e||) + Cn(g±)e||Cosh(L/eil)' 



where {q±) = y Cuiq±) / Cii {q±) / q± is the correlation length describing the decay of in-plane 
translational order in the z direction and cii(g_L) and C44{q±) are the non-local compressional 
and tilt moduli of the bulk flux array, respectively. The nonlocality of the elastic constants in 
the z direction is negligible at low fields compared to the in-plane variation. Finally, B2{q±) 
is the long wavelength compressional modulus of a 2d liquid of point vortices interacting via 
a l/r_L potential at large distances, so that B2{q±) ~ B'^/Airq^ as q± 0. As discussed in 
|TT| , translational correlations at the surface are controlled by the l/r±_ surface interaction 
only for small wavevectors, such that B2{q±) > CuL, or q± < q±^ = q±B2/ (Lcu). For 
> g_|_'^ the 2d surface structure factor is representative of that of a 2d cross-section of bulk 
and in the hydrodynamic model it is given by, 

nokBTqji'^ UokBT 

C44(g±)e|| Jc44(g±)cii(gx) 



The transfer function R{qj_,L) reduces then to 

R{q^, L) ^ [cosh(L/e||)]-' ^ 2e-^/^ii('?^), (6) 
where the second approximate equahty holds provided L » 61, i.e., if q± » q±* = 



C44/cii)/L. If the elastic constants are calculated from Ginzburg-Landau theory [jT5[, 
one finds q±^ ~ q±* ^ 1/L. On the other hand, there is evidence for a strong downward 
renormalization of the compressional modulus Cn from entropic effects at low fields [^0, 
as discussed below. In contrast the tilt modulus is expected to be accurately given by the 
Ginzburg-Landau theory. As a result, q±* >> 1/L. If the surface bulk modulus B2 is 
not renormalized, then q±^ = B"^ / {At^Lcu) ~ g_L* >> 1/L and surface effects control the 



surface translational order up to wavevectors of order g_L* ~ ~ 5/im ^, where we 



hand, we argue below that B2 may also be renormalized downward by entropic effects due 
to coupling of the surface tips to line wander in the bulk. In this case we expect B2 ~ cf^/g_L 
and g_L* ~ 1/L << q±*, so that the long-range surface interaction only controls surface 
translational correlations for wavevectors much smaller than those probed by the decoration 
experiments. 

The hydrodynamic model is very useful for describing the long-wavelength properties 
of the vortex array, and it allows us to incorporate the nonlocal effects of the intervortex 
interaction which are known to be important in flux crystals over much of the temperature- 
field phase diagram. An alternative more microscopic description can be obtained via the 
mapping of flux lines onto the world lines of two-dimensional bosons [Q. In the boson 
language the correlation length C,\\{q±) is determined by the Bogoliubov excitation spectrum 
e{q±) of a weakly interacting superfiuid, according to [^], 



where Vq = 0o/4^ = 47rA^eo is the energy scale of the bare intervortex interaction and ei is 
the single vortex tilt energy per unit length, with where A is the penetration length in the ab 
plane. Although ei ^ {M±/M^)eo\nK « eg for fields B » {M^/M-^ is the effective 
mass ratio), for 5 ~ 0/A^ (the regime relevant to decoration experiments) we have ei ~ eg 
because of magnetic couplings between the CUO2 planes [|14|. The result obtained from the 
analysis for the boson liquid agrees with the hydrodynamic result at small g_L, provided we 
make the identification C44 = noCi and cn = n^Vo. 

Translational order can be quite sensitive to point disorder, which is present in all exper- 
imental samples. The question of whether the flux patterns seen in decoration experiments 
are controlled by quenched disorder or by thermal fluctuations is at present open. Weak 
point disorder both in the bulk and at the surface of the sample can be incorporated in 
the hydrodynamic model as a random potential with short-range correlations coupled to 



have used the parameters of the Yao et al. experiment with L = 20fim [0]. On the other 
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the vortex density p. Bulk point disorder yields an additive Lorentzian squared correction 
to the thermal three-dimensional structure factor. Neglecting surface effects, the disorder 
contribution to the structure factor defined in Eq. |^ is given by []11 



SD{q±,L) = S2z,(gx)(l + Vei|)e-^/«", (8) 

where S2d{(1±) is the quenched-disorder contribution to the two-dimensional structure factor 
at one of the surfaces, 

and As is the correlator of the bulk random impurity potential. The total structure function 
is S{q±,L) = ST{q±,L) + SD{q±,L). The transmittance of translational order across the 
sample is governed by the same length scale C,\\ {q±) as in the thermal case independent of the 
strength of the quenched disorder. Weak surface disorder yields a contribution S2SD{q±) to 
the two-dimensional surface structure factor that vanishes as q±'^ at small wavevectors and 
can therefore be distinguished from the other contribuutions [|Tl] . 



Double-sided decoration experiments of the type carried out by Yao et al. can measure 
both the two-dimensional structure factor 5*2 (gj,) at one of the two surfaces, as well the 
correlations across the thickness of the sample described by S{q±, L). The transfer function 
R{q±,L) and then the "excitation spectrum" e{q±)/kBT are obtained from the ratio of 
these two correlation functions. Since the sample thickness is known, the slope of the 



Bogoliubov spectrum at small q± yields a measurement of the ratio ycu/cu. Yao et al. 
find C11/C44 ~ 1.5 X 10""^ in BSCCO single crystals at 12G, a value about four orders of 
magnitude smaller than predicted from the Ginzburg-Landau mean field theory [|15|. If one 
assumes that the main contribution to S2{q±) is the thermal one given by Eq. |^ and that 
the flux array falls out of equilibrium near the irreversibilty temperature Tj^r, i.e., Ty ^ Tj^r, 
one can also extract the geometric mean of the two elastic constants from the linear slope of 
5*2 (gj,) [|n|]. This gives C44 ~ 27G^ and Cn ~ 6 x lO^^G^ 0. The value of C44 is essentially 
equal to B'^/Ati and is consistent with C44 ~ ^0^0, provided one uses the value of A at the 



irreversibility line [jl^. The expression for the nonlocal tilt modulus obtained from the 



Ginzburg-Landau theory can be found in the Appendix of Ref. [O. If we use the expression 



for C44 that applies at low fields (ao < A, but q± >> ^ M^/ M^/ X so that nonlocal effects in 
the tilt modulus can be neglected), we obtain C44 ~ r2oeo/2, consistent with the experimental 
measurement. The experimental value for cn is about four order of magnitudes smaller than 
expected on the basis of Ginzburg-Landau theory, which neglects fluctuation effects. 

To obtain an approximate understanding of the strong downward renormalization of the 
compressional modulus (or equivalently of the strength of the intervortex interaction), we 
recall that the Bogoliubov results can be made quantitatively accurate for dilute superfluids, 
provided the bare interaction potential is replaced by an effective interaction or "t-matrix" 



deflned as the sum of an inflnite series of ladder diagrams |]T6| . In a two-dimensional super- 
fluid gas the renormalization of the q± = Q part of the intervortex interaction corresponds 
to the summation of a series in 1/ ln(l/noA^), where no is the boson density and A the range 
of the interaction, and leads to the replacement [^|[ 



^° 1 + [V,h/{kBm ln(l/noA2)/47r 

(10) 



ei ln(l/noA2)" 

The renormalized compressional modulus is then estimated as cf^ ~ n^Vn. Substituting the 
material parameters appropriate to the experiments of Ref. 0, we flnd cf^ ~ 2 x 10~^G^. 
The experiments may not be in the limit of extreme dilution {n^X^ « 1) required for 
the second line of Eq. so it is not surprising that this result is even lower than the 
experimental value. 

When nX^ > 0{1) one can qualitatively expect an analogous downward renormalization 



to arise from entropic contributions to the free energy from vortex-line braiding ||T8[. In this 
dense limit each vortex line spends a certain "time", i.e., length along the z axis, in the tube 
or "cage" of radius Oq provided by the repulsive interaction with its six neighbors. Flux 
lines wander within this cage until they escape to one of the approximately six neighboring 
cages. Collisions reduce the entropy of the interacting flux array relative to that of the 

8 



noninteracting system. Escape events, which yield flux-hne braiding, increase, however. 



the entropy, similar to the discussion of interstitial wandering in Ref. |]T8|. In the flux-line 
liquid, where escapes are frequent, the reduction in entropy due to collisions, or unsuccessful 
escapes, is a small correction. Each escape increases the entropy per vortex by fcslng, 
with q an effective coordination number describing the different directions in which a vortex 
can hop. The average distance Iz between hops among lattice sites is given by DqIz ~ a^, 
with Do = kBT/ii the vortex "diffusion constant" along z, or Iz ~ eia^/ksT In a 
sample of thickness L the total number of jumps is of order L/lz and the corresponding 
entropic contribution to the Gibbs free energy per unit volume of the vortex array is gent ~ 
~TL^^'^r •? = — -^^^p-n^ In q. The total Gibbs free energy per unit volume can be written 
as g{n) ^ -n^{H - H^i) + gint{n) + gent{n), where gint ~ eon^A^ is the contribution 
from intervortex interactions. Upon expanding about the minimum density to otain cn = 
{(Pg/dn'^)\n=no, "we see that the entropic contribution partially cancels the large contribution 
from interactions, consistent with experimental observations. The entropic and energetic 



contributions to Cn are comparable when T ^ veoIiA, which is comparable to the melting 
temperature of the Abrikosov flux lattice in this field regime |p.8|| . 

An analogous mechanism could lead to a strong downward renormalization of the surface 
interaction. This is because the flux tips at the sample surface are not true point vortices, 
but are connected to the flux lines in the bulk. Braiding effects of the type described above 



within a surface layer of thickness ^zic[±) ~ y cu/cii/q± will increase the surface entropy 
of the fiux-tips, yielding a contribution g^^^ ^ —kBTn{S,z/lz) Ing ^ — -^^|=^n^/^ In g to the 
free energy per unit area, where we have used C44 ~ nli and cn ~ n^Vo. The corresponding 
free energy from surface interaction among the flux tips is gint ~ ""-^i^^- Provided ei ~ eo, 
as is appropriate for this low field regime, the two contributions to the energy (and hence 
to B2{q±)) are again comparable near the melting temperature. 

The finite range of the intervortex interaction can be incorporated in the derivation of the 
Bogoliubov spectrum, which is then given by Eq. |^, with Vq = P,pi? 



The second term in Eq. |^ is unchanged since it represents the "kinetic energy" contribution 



to the spectrum, which is unrenormahzed due to Gahleian invariance of the equivalent 
boson problem. The full wavevector-dependent interaction V{q±) is again renormalized by 
resumming an infinite series of ladder diagrams, which leads to an integral equation for 



the effective t-matrix at finite wavevector W^. Upon neglecting corrections logarithmic in 



the wavevector, we find that the screening length A is not renormalized and one obtains 
VR{q±) = Vr/{1 + g±^A^). The Bogoliubov spectrum can then be rewritten in a suggestive 
form that interpolates between the boson result and the hydrodynamic description as 



" ksT V C44 ^ 2ei ^ 



where we have identified the renormalized local compressional modulus as cf^(g_L) = 
nQVfi{q±). The nonlocality of the tilt modulus is not important for the low fields of in- 
terest here |jl5|, as confirmed by the experimental finding that C44 ~ ^o^i- The renormalized 



Bogoliubov spectrum given in Eq. |TT] is shown in Fig. 1. 

The Bogoliubov spectrum is not expected to be quantitatively accurate for dense super- 
fluids. In this regime the theory can be improved following Feynman and approximating 
this spectrum by 

ksT 2nQliS2{q±y 

where 5*2 (gj,) is the structure factor of a two-dimensional cross section of a dense vortex 
liquid. Its thermal contribution in hydrodynamic theory is given in Eq. |[ For more realistic 
functions 5*2 (gj,), this formula leads to a "roton" minimum in the excitation spectrum at 
q± ^ ksz = v^47m^, at approximately the position of the first maximum in S2{q±). This 
"roton" minimum has been observed in the experiments |^. 

III. TRANSLATIONAL AND ROTATIONAL REGISTRY 

A source of uncertainty arises in the experiments from the difficulty in matching the (x, y) 
locations of vortices being imaged on the two sides of the sample. This positional uncertainty 
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can be decreased when the sample contains locahzed defects that run all the way across the 
sample, such as grain boundaries, since these can provide a common reference frame on the 
two sides 0. Correlations on different length scales are in general affected differently by this 
mismatch. To quantify this effect for a given sample we define the rms density fluctuations 
arising from a translational mismatch d and an orientational mismatch (p averaged over the 
area A of the sample, as 

72 

A{d, (t),L)= j [5n{n^ • + d, L) - 5n{v^, Q)f >, (13) 

where TZ^ is a two-dimensional rotation matrix and we neglect effects due to quenched 
random disorder. We expect A((i, 0, L) to be a minimum when the patterns on the two 
sides are defined relative to (x, y) coordinate systems with a common origin and orientation. 
Minimizing A with respect to d and using experimental data could be used to bring 
these coordinate systems into registry even in the absence of identifying features such as 
grain boundaries which penetrate across the entire crystal. Upon introducing the Fourier 
components of the density, Eq. |1^ can be rewritten using the single-pole approximation in 
g-space for the structure factor. 



as. 



A{d, 0, L) = / c?g^g^3^_ 

27rC44 Jo ^r(Q\ . 

X 



i_j„(,,,)Zl(|?^5£!^e--.<-V'=.-], (15) 

RkBz sm(0/2) J 
where JqI^^) ^^nd Ji(x) are Bessel functions and R denotes the linear dimensions of the sam- 
ple in the ah plane, with A = vri?^. The integral on the right hand side of Eq. |l^ diverges 



and is cutoff by a circular Brillouin zone, ksz = V^™o- We have evaluated the dimen- 
sionless quantity A{d,(f),L) = [A((i, 0, L) — A(0, 0, L)]/A(0, 0, L) using the hydrodynamic 



approximation eR{q±)/kBT = q±\/cii{q±)/c44. The function A{d,(j),L) is shown in Figs. 
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2 as a function of both the angle and the translation d for a few values of the sample 
thickness. The function A((i, 0, L) has a parabolic minimum at ci = 0, = 0, according to 



A(d, 0, L) ^ A(0, 0, L)^ 1 + ^a(L) [^^4^/2 + sin2(0/2) 



(16) 



where the dimensionless curvature a{L) is given by 



a{L) 



L*kBz- 



BZ 



-L*kf 



1-e 



(17) 



and L* = L^cii{q± = 0)/c44. At low density, corresponding to Xksz << 1, the behavior 
of the curvature is controlled by L*kBz, with a ~ \/{L*kBz) for L*kBz « 1 and a ~ 
l/{L*kBzY L*kBz >> 1- At high density, corresponding to XkBz >> 1; the relevant 
length scale is L^A and a ~ X/L* for L* « X and a ~ l/{L*kBz)^ for L* >> A. 



IV. TRANSMISSION OF ORIENTATIONAL ORDER 

From the analysis of double-sided decorations one can also study the propagation of orien- 
tational order across the sample. Orientational order is much less sensitive to point pinning 
PCI|| . It is measured by correlations in the bond-orientational order parameter ^p6{r) = e^^^^^\ 
where ^(r) is the bond- angle field. The corresponding angular correlation across the sample 
thickness is 

^ exp[-18 < [9{r^, L) - 0(0, 0)]^ >]. (18) 

The decay of bond-orientational order in a hexatic flux liquid was discussed in Ref. |Q] 
in the hydrodynamic limit. In a bulk sample, ignoring boundary conditions and surface 
effects, the in-plane Fourier transform of the thermal part of the correlation function of the 
bond-orientational order parameter was found to be given by 

GH{q±,L) =< It/^sT > 
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^5qi,o + GH2{q±)e 



(19) 



with 



9kBT g/cijT 1 

GH2(g±) = -r^^H{q±) = I , — ■ (20) 



Here ^h{(1±) = y K\/ K\/ q_[_ is the correlation length governing the transmittance of hexatic 
order across an hexatic flux-line liquid and K\ and-ft';4 are the hexatic stiffnesses. 

In a superconducting slab of flnite thickness L we use free boundary conditions on the 
bond-angle fleld at the surface to flnd that Eq. |19| is replaced by, 

GuiqL, L) =< l^el' > [^5qx,o + GH2iq±)RH{q±, L)] , (21) 
with GH2{q±) = '^^Hiq±)coth{L/^H) and Rniqi^L) = [cosh(L/^^,)]-i. 

A 

In a flux lattice with long-range crystalline order the hexatic order parameter is not 
independent, but is simply related to the curl of the elastic diplacement fleld, 6* = |z • 
(V X -u). The correlation function of the bond-angle fleld in an inflnite sample then follows 
immediately, with the result, 

G^(gx, L) =< l^eP > [^^^,,0 + G|,2(g±)e~^/^-^'"^] , (22) 

where ^H{q±) = ^J^/q± is the correlation length governing transmittance of hexatic order 
across a flux lattice, and 



Gf,M = TT^iffrT = T^-Jx. (23) 



4c44^i5'(g±) 4^C66C44' 

The corresponding expressions in a flnite-thickness sample with free boundary condition on 
the bond-angle fleld are modifled with the same flnite-size functions of L/C,H{qj_) as in the 
case of the hexatic flux liquid. This result shows that by measuring the correlation of bond 
order across the sample, as well as 6*^2 (?J-) at one of the surfaces, one can infer the value of 



the tilt and shear moduli. In addition, Eqs. ^ and ^ show that surface bond-orientational 
order decays as l/g^ in a hexatic liquid, but grows as q± in a lattice, providing a further 
mean to distinguish hexatic and crystalline order in the vortex array. 

Bond order decays exponentially in flux liquids which are isotropic in a plane perpen- 
dicular to the fleld direction. The results in this case are similar to Eq. |19|, except that the 
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delta-function term is absent and lim^j^-^o Ch{(1±) is finite and equal to the hexatic correlation 
length along the z direction. 
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Figure Captions 

Fig. 1. The renormalized Bogoliubov spectrum given in Eq. ( ]Tl| ) as a function of wavevector 
for uqX'^ = 0.2,0.5,1. 

Fig. 2. The spatially averaged mismatch function A{d, (p, L) is shown (a) at = as a 
function of d and (b) at c? = as a function of 0, for three values of L. Note that the 
sample thickness only enters in the dimensionless combination LkBz\J cii{q± = 0)/cm. We 
have used B = 12G, cii(g_L = 0)/cu = 1.5 x 10"'^ and R = 0.2mm. 
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